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Abstract 

In this paper we will discuss an algorithm for extracting skewed parton 
distributions from experiment as well as the relevant process and experimental 
observable suitable for the extraction procedure. 
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I. INTRODUCTION 

Skewed parton distributions^ which appear in exclusive, hard diffractive processes like 
deeply virtual Compton scattering (DVCS)0 or vector meson production with a rapidity gap, 
to name just a few, have attracted a lot of theoretical and experimental interest over the 
last few years as a hot bed for interesting new QCD physics 0-^,0. The list of references 
is probably far from being complete and thus we apologize beforehand to everybody not 
mentioned. 



This is the unified terminology since the Regensburg conference of '98, finally eradicating the 
many terms like non-diagonal, off-diagonal, non-forward and off-forward which have populated the 
literature on this subject over the last few years. However recent publications have, alas, again 
fallen back upon the old terminology! 
2 First discussed in Ref. (l|. 
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FIG. 1. The lowest order handbag contribution to DVCS with Q 2 = —q 2 and q' 2 = 0. 

The basic concept of skewed parton distributions is illustrated in Fig. [I] with the lowest 
order graph of DVCS in which a quark of momentum fraction x\ leaves the proton and is 
returned to the proton with momentum fraction x%. The two momentum fractions not being 
equal is due to the fact that an on-shell photon is produced which necessitates a change in 
the + momentum in going from the virtual space-like photon with + momentum usually 
taken to be — xp + , where p + is the appropriate light cone momentum of the proton and x is 
the usual Bjorken x, to basically zero + momentum of the real photon. This sets 22 = X\ — x 
and thus the skewedness parameter to x. (see Q for more details on the kinematics.) 

Thus one has a nonzero momentum transfer onto the proton and the parton distributions 
which enter the process are non longer the regular parton distributions of inclusive reactions 
since the matrix element of the appropriate quark and gluon operators is now taken between 
states of unequal momentum rather than equal momentum as in the inclusive case (see for 
example ||). These parton distributions still obey DGLAP-type evolution equations but of 
a generalized form (see for example Radyushkin's references in 0). 

The above mentioned kinematical situation is not the only one possible. One can also 
have the situation where x<i becomes negative. In this case not a quark is returned to the pro- 
ton but rather an anti-quark is emitted. In this situation one does not deal with parton distri- 
butions any more but rather distributional amplitudes obeying now Efremov-Radyushkin- 
Brodsky-LePage (ERBL) type evolution equations (again see, for example, Radyushkin's 
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references in @). Furthermore, both momentum fractions could be negative in which case 
one is dealing with anti-quark distributions which again obey DGLAP-type evolution equa- 
tions. 

After having answered the question how these skewed parton distributions arise, the next 
question is which of the exclusive, hard diffractive processes is most suitable for extracting 
these skewed parton distributions and how can this be achieved. This question will be 
answered in the following sections, where we discuss the most promising process and the 
appropriate experimental observable in Sec. fH, in Sec. HJ we will explain the algorithm and 



the problems associated with it and finally in Sec. [TV] we will give an outlook on further 
research in this area. 



II. APPROPRIATE PROCESS AND EXPERIMENTAL OBSERVABLE 

The most desirable process for extracting skewed parton distributions is the one with 
the least theoretical uncertainty, the least singular Q 2 behavior so as to be accessible over a 
broad range of Q 2 and with a proven factorization formula. The last requirement is actually 
the most important one since without a factorization theorem one has no reliable theoretical 
basis for extracting parton distributions. 

The process which fulfills all the above criteria is DVCS since it is least suppressed in Q 2 
of all known exclusive, hard diffraction processes, in fact it is only down by an additional 
factor of Q 2 in the differential cross section as compared to DIS0, the theoretical uncertainty 
is minimal since we are dealing with an elementary particle in the final state as compared 
to, for example, vector meson production where one also has to deal with the vector meson 
wavefunction in the factorization formula as an additional uncertainty and there exists a 
proven factorization formula ||. 



3 Compare this to the 1/Q 8 behavior of vector meson production, di-muon production or di-jet 
production. 



Furthermore it was shown in Ref. |4j] that there will be sufficient DVCS events at HERA 
as compared to DIS, albeit only at small x between 10~ 4 — 10 -2 , to allow an analysis with 
enough statistics. 

The experimental observable which allows direct access to the skewed parton distribu- 
tions is the azimuthal angle asymmetry A of the combined DVCS and Bethe-Heitler(BH)Q 
differential cross section, where the azimuthal angle is between the final state proton - 7* 
plane and the electron scattering plane. A is defined as [[§]: 

^ f-i%d(f) da DV cs+BH — f^ji" 2 d4> da DV cs+BH 
J* d(f> daovcs+BH 

In other words one counts the events where electron and photon are in the same hemi- 
sphere of the detector and subtracts the number of events where they are in opposite hemi- 
spheres and normalizes this expression with the total number of events. 

The reason why this asymmetry is not is due to the interference term between BH and 
DVCS which is proportional not only to cos(0), as compared to the pure DVCS and BH 
differential cross sections which are constant in 0, but also to the real part of the DVCS 
amplitude. The factorized expression for the real part of the amplitude takes the following 
form 0: 

Re T(x, Q 2 ) = f 1 ^-Re Q(x/y, Q 2 )f t (y, x, Q 2 ). (2) 
J-i+x y 

Re Ci is the real part of the hard scattering coefficient and fi are the skewed parton dis- 
tributions. The sum over the parton index % is implied and y is defined to be the parent 
momentum fraction in the parton distribution. As mentioned above, one is mainly restricted 
to the small- 2 region where gluons dominate and thus i in Eq. (JJ) will be only g to a very 
good accuracy. Note that since the parton distributions are purely real, the real part of the 
amplitude in its factorized form has to contain the real part of the hard scattering coefficient. 



4 In the Bethe-Heitler process the incoming electron exchanges a coulomb photon with the proton 
and radiates off a real photon, either before or after the interaction with the proton. 
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Thus Eq. (HD contains only measurable or directly computable quantities except Eq. (0) 
in the interference part of the differential cross section for real photon production which is 
isolated in Eq. ([]]). Therefore, one would now be able to extract the skewed parton distribu- 
tions from experimental information on A, the directly computable part of the interference 
term and the knowledge about the hard scattering coefficient if one could deconvolute Eq. 

(!)• 

As we will see in the next section this direct deconvolution is not possible, however there 
is a way around the deconvolution problem. 



III. ALGORITHMS FOR EXTRACTING SKEWED PARTON DISTRIBUTIONS 

A. The Deconvolution Problem in DIS 

The deconvolution problem in inclusive DIS presents itself in a similar way as in Eq. 
(0). For the structure function F 2 (x,Q 2 ), for example, one has the following factorization 
equation (in a general form): 

F 2 (x,Q 2 ) = f'^-Qix/y^^f^Q 2 ), (3) 
Jx y 

where one has the same situation as in Eq. (0) except, of course, that the hard scattering 
coefficient Cj and the parton distributions fi are now different from the DVCS case. Also 
notice that the parton distributions depends now only on y rather than y and x. In this 
case one can now easily deconvolute Eq. (|3]) by taking moments in x via is an 

easy exercise to show that the convolution integral turns into a product in moment space: 

F 2 (N,Q 2 ) = C l (N,Q 2 )f l (N,Q 2 ). (4) 

Thus, after having calculated the hard scattering coefficient to the appropriate order 
and having measured F 2 such that the moment integral can be taken numerically, one can 
directly extract the parton distribution. What remains to be done is to perform the inverse 
Mellin transform to obtain the parton distribution in terms of x and Q 2 . Of course, we have 
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simplified the actual procedure and the inverse Mellin transform is also not easy to perform 
but this example serves more as a pedagogical exercise to illustrate the basic concept of 
deconvolution and extraction of parton distributions. 

In the case of interest to us, however, life is not that "simple", since the skewed parton 
distributions depend on two rather than one variable. Furthermore, the hard scattering 
coefficient depends on the same variables as the parton distribution. This makes the decon- 
volution of Eq. @, at least to the best knowledge of the author, impossible because both the 
hard scattering coefficient and the parton distribution have two rather than one variable in 
common, one of which is even fixed, thus one does not have enough information to perform 
a deconvolution. 

This seems like an intractable problem but there is a way out. For the purpose of as 
simple a presentation as possible the following discussion will only be done in LO but the 
same principles also apply in NLO. However, the precision of the data in the foreseeable 
future, will be such that a leading order analysis will be sufficient. The following two 
discussions rest heavily on the methods in Ref. 

B. The First Principle Extraction Algorithms 

The basic idea of this algorithm is to expand the parton distributions in terms of orthog- 
onal polynomials to reduce the unknown quantities in the factorization formula for the real 
part of the DVCS amplitude to a number of unknown coefficients which can be obtained 
through an inversion of a known matrix and DVCS data on the real part of the amplitude. 

As is well known, one can expand parton distributions or any smooth function for that 
matter, with respect to a complete set of orthogonal polynomials Pj ap \t), where t is used 
here to shorten the notation. The orthogonality of the polynomials of our choice needs 
to be on the interval — 1 < t < 1 with t = which translates to an interval in y of 
— l+a;<y<las found as the upper and lower bounds of the convolution integral in Eq. 
(0). One can then write the following expansion: 
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/M(t ,,, Q2) = _l_g^Mp r(t) M r( ,,Q^) (5) 

with w(t\acp) and rij(ap) being weight and normalization factors determined by the choice 
of the orthogonal polynomial used. The labels q, g for quarks and gluons are necessary since 
the j label will be different for quarks and gluons. «p is a label which depends on the 
orthogonal polynomials usedF]. Mj ,9 (x,Q 2 ) is given by: 

oo 

Mf 9 {x, Q 2 ) = J2 Ef k 9 {v- a P \x)f k > 9 (x, Q 2 ), (6) 

fc=0 

where 

f k 9 (x,Q 2 ) = j:x l - k B? k *fr(x,Q 2 ). (7) 

1=0 

Bfjf is an operator transformation matrix which fixes the NLO corrections to the eigenfunc- 
tions of the kernels. The explicit form of the transformation matrix B q k can be found in 
Eq. (35) of the second article of Ref. || for example. The explicit form is not important, 
however, for neither this discussion nor the conclusion of this paper since we are dealing 
only with a LO analysis in which case the transformation matrix is just the identity matrix. 
The general from was just included for completeness sake. 

The moments fi' 9 {x, Q 2 ) of the parton distributions in Eq. (^) generally evolve according 

to 

fr(x, Q 2 ) = Ki k (a s (Q 2 ), a s {Ql))tf' 9 {^ Ql) (8) 

where the evolution operator is a matrix (i, k equals either q or g) of functions in the singlet 
case (and just a function in the non-singlet case) taking account of quark and gluon mixing 
and depending on the order in the strong coupling constant. 

Striving for simplicity, the above expansion is most simple in LO in the basis of Gegen- 
bauer polynomials, since they are the eigenf unctions of the evolution kernels at LO. Thus 



5 ap = a, (3, in other words two labels, if Jacobi polynomials are used or ap = fi— 1/2 if Gegenbauer 
polynomials are used. 
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we will use these polynomials from now on in our formulas. Thus, the Gegenbauer moments 
of the initial parton distributions at Ql in Eq. (g) are defined the following way 

mx,Qi) = J 1 _ i dt (^) l cf /2 (^)nt,x,Qi) 

f*(x,Qt) = \[dt {^r^Ctl (^) f 9 (t,x,Ql). (9) 

Turning now to the expansion coefficients in Eq. (H) . The upper limit of the sum in Eq. 
(0) is given by the constraint ^-functions Q present in the expansion coefficients, which are 
defined, in terms of Gegenbauer polynomials, by 

i» (-i)^r(/z + ffi 
i» r> + fc)r(i + ^) 

where z/ = /x = 3/2 for quarks and z/ = /x = 5/2 for gluons. 
The general expansion in Eq. (J|) then reduces to 

nt, x,Q 2 ) = ^tf w^ E umx)cf\t)ft{ X , g 2 ) (ii) 

A X j=0k=0 Iv j\ 6 / Z ) 

f 9 (t, ^,Q 2 )~EE ^f^^-i(5/2|x)Cf ^t)/!^^, Q 2 ), (12) 
z x j=0k=l Iy j\°/ Z ) 

with iu(t|i/) = (t(l - t)) 1 '" 1 / 2 , xVj(i/) = 2- 2 ^+ i r2 r iy ) ^+f ) and the are Gegenbauer 
polynomials. 

As we are in LO, multiplicatively renormalizable moments evolve with the following 
explicit evolution operator: 

Kf (« S (Q 2 ), a s (Ql)) = Texp yjf(a s (r))) (13) 

where T orders the matrices of LO anomalous dimensions along the integration path. Note 
that there is a slight difference in the anomalous dimensions in the skewed case to the 
anomalous dimensions in the non-skewed, i.e. inclusive, case due to the particular definition 



% k = l, if k<j; 0, iij<k 



i + (-i 



\3-k 
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of the conformal operators used in the definition of the parton distributions^] 7J 9 = |^/| ff '* nc ■ 
and 7 f = tr'"''- 

Now we have all the ingredients to proceed. Inserting Eq. fll2"|) in Eq. (0) one obtains for 
small x, where one is justified to neglect the quark contribution: 

00 00 

Re T(x, Q 2 ) = 2 £ £ K 9 k 9 _ x (a s (Q 2 ) , a.{(&))fk-i(*, Ql)E%-^/2\x) 

j=0 fc=l 

7-i 2t + x Nj{5/2) 9 \2 x J J ~ lW ' V ; 

where we chose the factorization scale to be equal to the renormalization scale , which we 
chose to be equal to Q 2 . As one can see the integral in the sum is now only over known 
functions and will yield, for fixed x, a function of j as will also the expansion coefficients for 
fixed x. The evolution operator can also be evaluated and will yield for fixed Q 2 also just 
a function of j, which leaves the coefficients fl-i{ x iQl) as the only unknowns, albeit an 
infinite number of them. Since the lefthand side will be known from experiment for fixed x 
and Q 2 , we are still in the unfortunate situation that a number is determined by the sum over 
an infinite number of coefficients labeled by j. Thus, if one had measured the real part of the 
DVCS amplitude through the asymmetry A at fixed x | and at an infinite number of Q 2 , one 
would have an infinite dimensional column vector on the lefthand side namely the real part 
of the amplitude at fixed x but at an infinite number of Q 2 and on the right hand side one 
would have a square matrixP] times another column vector of coefficients of which the length 



7 This is true in LO, in NLO, however, the anomalous dimensions obtain, besides the NLO 
jj's of the inclusive case, additional anomalous dimensions due to non-diagonal elements in the 
renormalization matrix of the conformal operators entering the skewed parton distributions (see 
Ref. §). 

8 This fixes the undetermined coefficients up to the j index. 

9 The number of Q values determines the column dimension and the number of the index j deter- 
mines the row dimension. The matrix is square since we can choose the number of Q values to be 
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is determined by the number of j. Since all the entries in the matrix are real and positive 
definite^ it can be inverted, using the well known linear algebra theorems on inversion of 
infinite dimensional square matrices, provided that there are no zero eigenvalues in other 
words no physical zero modes in the problem which would imply that the real part of the 
DVCS amplitude would have to be zero which is, of course, never the caseQ. After having 
found the inverse, we can directly compute the moments of our initial parton distributions^ 
which are needed to reconstruct the skewed gluon distribution at small x from Eq. (|5|). 

The drawback of the above procedure is that this process has to be repeated anew for 
each x. Nothing, however prevents us from doing so, in principle. Even for a finite number 
of j's and Q's, the task seems formidable, however, this is not as problematic as it seems, 
since experiment will only render information for small x, at least in the beginning, and not 
the whole range of x, thus one does not need an infinite number of coefficients and thus an 
infinite number of Q for each x to get a good approximation. Unfortunately aj max of 50 — 100 
will be necessary, 0. therefore, if the lefthand side is known for each x at 50 values of Q 2 , 
Eq. reduces to a system of 50 equations with 50 unknowns for j max = 50. This system 
can readily be solved as explained above. Experimentally speaking, of course, this procedure 



equal to the number of j values! 

10 The evolution operator will, of course, always yield a positive number, the integrals in the 
sum, are integrals over positive definite functions in the integration interval and the expansion 
coefficients are also positive definite as can be seen from Eq. fllPp, 

11 The DGLAP part of the amplitude will be zero at x = 1 but the contribution from the ERBL 
region will not be! 

12 Note that the same moments of the initial parton distribution will appear for different values of 
j, since the sum over k runs up to j, for fixed Q, such that each unknown moment is just multiplied 
by a number determined from known functions in Eq. (|l4|). 

13 The author's thanks go to Andrei Belitsky for pointing this out. 
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is not feasible, though theoretically very attractive, since one will never be able to measure 
the any experimental observable for fixed x at 50 different values of Q 2 . Nevertheless, there 
may be ways using constraints on SPD's to reduce this number of 50 — 100 polynomials as 
can be done in the forward case, however this has to be further explored. 

Notwithstanding the above, let us give a toy example of the above extraction algorithm 
0. Take x discrete and fix Q 2 then one can write a factorized expression for a cross section: 

°a = Y. H r,afr,a- (15) 

j 

The index j corresponds to the parton fractional momentum, and a to the x variable. 
Obviously, it is not possible to obtain fj. a from a a . If one now puts in an index for Q, the 
parton densities will now be f(Q)j ;a , and the solution of the evolution equation has the form 

f(Q)r,a = T, U (Q)^af(Qo) k] a. (16) 
k 

Here f(Qo)k;a is the initial parton density at the value Q = Qo which is left implicit. The 
cross section as a function of Q takes now the form 

a Q;a — ^ Hj;aU(Q)j,k;af(Qa)k;a 
j,k 

= J2 A (Q)j,kf(Qo)k;a, (17) 
k 

for a suitable matrix A. The Q dependence of the hard scattering function H can be ignored 
for our present purpose. 

As a next step, take enough values of Q such that the matrix is square. The most trivial 
example is to have two values of Q: the initial value and one other: 

fi-j — f(Qo)j 

/ 2; ;=EWW(QoV (18) 
k 



The author would like to thank John Collins for suggesting such an example to clarify the 
problem at hand. 
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One can take U to be triangular, as is appropriate for DGLAP evolution 

„=(; ;), 

the hard scattering cross section to be 

#=(1,1) (20) 



and the parton distributions to be a two dimensional column vector: 

f(Q 
AQo)2 



, f(Qo)i , 

fo = . (21) 



This then yields 

<r(Q) \ _n 2\ //(Qo)i 

,(7(Qo)/ \1 1/ U(Qo)2 



(22) 

Clearly one has an invertible matrix in Eq. (|22] ) and can thus compute f(Qo)i and f(Qo)2- 



C. The Practical Extraction Algorithm 

A practical way out of the polynomial predicament is by making a simple minded ansatz 
for the skewed gluon distribution, since we are still at small x, in the different ER-BL and 
DGLAP regions of the convolution integral Eq. (§). An example of such an Ansatz could 
be Aqz~ Ai {1 — 2)3 for the DGLAP region where z is now just a dummy variable. If one 
inserts this Ansatz in Eq. (fj) and can fit the coefficients to the data of the real part of the 
DVCS amplitude for fixed x and Q 2 . One can then repeat this procedure for different values 
of Q 2 and then interpolate between the different coefficients to obtain a functional form of 
the coefficients in Q 2 . Alternatively, after having extracted the values of the coefficients for 
different values of x at the same Q 2 , use an evolution programm with the ansatz and the 
fitted coefficients as input and check whether one can reproduce the data for the real part 
at higher Q 2 , thus checking the viability of the model ansatz. 

To obtain an ansatz fulfilling the various constraints for SPD's (see Ji's and Radyushkin's 
references in PJ|), one should start from the double distributions (DD) (see Redyushkin's 
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references in which yield the skewed gluon distribution in the various regions of the 

convolution integral 



Due to the fact that there are no anti-gluons, the above formula is enough to cover the whole 
region of interest — 1 + x < yleql. What remains is to choose an appropriate model ansatz 
for G, for example, 



with f(zi) being taken from a diagonal parametrization with its coefficients now being left 
as variants in the skewed case and the normalization condition h(zi) = Jq~ z1 dzh(zi, z) such 
that, in the diagonal limit, the DD just gives the diagonal distribution. The choice for 
h(z±, z) is a matter of taste but should be kept as simple as possible. The drawback of this 
algorithm as compared to the previous one is that it is model dependent and thus not a first 
principle methods, which theoretically speaking, is not as satisfying but from the practical 
side this method is much simpler and thus experimentally much more feasible. 

Thus, one has solved the problem of extracting the parton distributions from the factor- 
ization equation, at least for small-x. The remaining problem is an experimental one. 



After having showed, that the extraction of skewed parton distributions from DVCS 
experiments is both principally and practically possible given the high enough statistics data 
on the asymmetry, one should now get a more accurate model description of the asymmetry. 
This will be done elsewhere. 
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